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ABSTRACT 

Experimental and numerical evidence is reviewed for the existence of a Stewartson 
layer in spherical Couette flow at small Ekman and Rossby numbers (E 10~ 3 , 
Ro ;$ 10 -2 ), the relevant hydrodynamic regime in the superfluid outer core of a 
neutron star. Numerical simulations of a superfluid Stewartson layer are pre- 
sented for the first time, showing how the layer is disrupted by nonaxisymmetric 
instabilities. The unstable ranges of E and Ro are compared with estimates of 
these quantities in radio pulsars that exhibit glitches. It is found that glitching 
pulsars lie on the stable side of the instability boundary, allowing differential 
rotation to build up before a glitch. 

Subject headings: dense matter — hydrodynamics — stars: interior — stars: 
neutron — stars: rotation 



1. Introduction 

Meridional circulation, driven by Ekman pumping, occurs routinely in the atmospheres, 
oceans, and fluid interiors of rapidly rotating astrophysical objects. Indeed, it is a generic 
feature of Navier-St okes flow in any spheri cal Couette geometry (i.e., a differentially rotating, 



spherical shell); see lJunk fc Egbersl (120001 ) for a review. As the Ekman number E decreases, 



and the differential ro tation increases, spherica l Couette flow becomes nonaxisymmetric and 



eventually turbulent (jNakabayashi et al.ll2002l ). In the limit of rapid overall rotation, a de- 



tached shear layer, known as the Stewartson layer, forms along the tangent cylinder to the 
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inner sphere (jStewartsonl 119571 . Il966l ; iBussd Il968l ) . It can be disrupted by nonaxisymmet 



ric instabilities, 
thereby excited ( 


Numerical simulations indicate that 


Hollerbach 


1994; 


Dormv et al. 


1998 


2005; 


Hollerbach et al. 


2006 


. 



The possible existence of an unstable Stewartson layer in a differentially rotating neu- 
tron star has important astrophysical consequences. This is true especially if the inner core 
rotates faster than the rest of the star, like in the Earth, a real possibility in models where the 
inner core makes a transition t o a crystalline color-superconducting phase (lAlford fc Reddy 
20031 ; lAlford et al.l 120051 . 120081 ). Recently, the importance of the global flow pattern inside 
a neutron star to the phenomenon of puls ar glitches has been highlighted by simu lations 
on the vortex (IWarszawski &: Melatod 120081 ) and hydrodynamic (IPeralta et al.ll2005l ) levels. 
Observational data sugg est that glitches resu lt from scale-invariant vortex avalanches driven 
by differential rotation (IMelatos et aJj 120081 ) . If the meridional circulation is fast enough, 
a vortex tangle is alternately created and destroyed, producing impulsive and oscil latory 
torque variations (IPeralta et al.ll2006l ; lAndersson et al.ll2007l ; IMelatos fc Peraltal 120071 ). The 
presence of a Stewartson layer modifies these conclusions and those of other studies (e.g., of 
stellar oscillations), w here a multicqmponent superfluid is perturbed starting from a nontriv- 
ial eq uilibrium state (jPeralta et al. 2006 . 2008 ; Glampedakis et al. 2009 ; Passamonti et al. 
2009h . 



To date, no studies have been published of the formation and stability of a Stewartson 
layer in superfluid spherical Couette flow. In this Letter, we present the first numerical 
simulation of such a system as an idealized model of the superfluid outer core of a neutron 
star. We calculate stability curves for a range of unstable nonaxisymmetric modes and 
compare the conditions for instability with observational glitch data, finding an upper limit 
on the velocity shear and hence the glitch sizes observed. The paper is organized as follows. 
In §5J we briefly review the analytic and numerical theory of Stewartson layers in viscous 
fluids, before calculating the structure of a steady Stewartson layer in a 1 So-paired neutron 
superfluid. In §3], we study the stability of the layer to nonaxisymmetric perturbations as a 
function of Ekman number E and Rossby number Ro. In §H we compare the stability basin 
in the E-Ro plane with available glitch data. The astrophysical implications are discussed 
in m 



2. Stewartson layers in neutron stars 



Consider a viscous fluid flowing inside a differentially rotating spherical container, with 
inner radius (angular frequency) Ri (f^i), outer radius (angular frequency) R 2 (^2), Rossby 
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number Ro = ^i)/^2 *C 1, and Ekman number E = v n j (i?2~ -Ri) 2 ^2, where v n denotes 
the kinematic viscosity. In a frame corotating with the outer sphere, the fluid outside the 
cylinder tangential to the in ner sphere is at r est, while the fluid inside the tangent cylinder 
moves in a columnar fashion (|Proudmanlll956l ) . Fluid is expelled (sucked in) by Ekman layers 
at r = i? 2 (Ri), while a triple-deck Stewartson layer buffers the jump in angular velocity 
across the tangent cylinder. It consists of an inner layer of thickness E 1 ^ 3 sa ndwiched betwee n 
layers of thicknesses E 2 ^ (E 1 ^) just inside (outside) the tangent cylinder (jStewartsonlll966l ) . 
The Ekman layers scale as E 1 / 2 , except near the equator of the inner sphere, where they 
scale as E 2 / 5 . 

A superfluid Stewartson layer in a spherical shell exhibits a similar structure. Figure 
QJi graphs the angular velocity in the rotating frame as a function of the cylindrical radius 
s = rsin#, for E = 1 x 10~ 3 (upper curve), 1 x 1CT 4 (middle curve), and 2 x 1CT 5 (lower 
curve), with Ro = 10~ 4 . The layer starts at cylindrical radius s ~ 1.8 and its thickness 
decreases with decreasing E, extending out to s ~ 2.7 for E = 1 x 10~ 3 and s ~ 2.1 
for E = 2 x 1 0~ 5 . In viscous fl ows, the thickness of the layer changes by less than 1 % for 
r$ Ro £j 0.5 (]Hollerbach!l2003l ) : similar behaviour is observed here. The inner Ekman layers 
are thicker at the equator. 

Figure [lb displays meridional streamlines for E — 1 x 10~ 3 (left), 1 x 10~ 4 (center) and 
2 x 10~ 5 (right), with Ro = 10~ 4 . The Stewartson layer is visible along the tangent cylinder, 
narrowing from left to right. The streamlines are drawn in the rotating frame of the outer 
sphere; the blank region to the right of the tangent cylinder indicates that the fluid there is 
at rest. The characteristic meridional speed in the layer scales as E 0A (R2 — Ri)Q2- 

To obtain the results in Figure [lj we solve the two-component Hall-Vinen-Bekarevich- 
Khalatnikov (HVBK) equations for a superfluid inside a spherical differentially rotating shell, 
with ffi / Ro = 0.67, using a pseudospectral collocation and time-split method (IPeralta et al. 



20051 . I2008T ) . The details of the calculation will be set out in a longer paper. Bound- 



ary conditions assume the presence of an inner core or a tra nsition between a 1 S'o and 
3 P2 superfluid ( Yakovlev et al. 1999 ; Mastrano Sz Melatosl 2005 ). We ignore vortex pinning 
and prot on- neutron entrainment for s i mplicity, although recent work shows it to be im- 
portant (ISedrakian fc Sedrakianl Il995l ; lAndersson fe Corner! l200ll ). We adopt no-slip and 
no-penetration boundary conditions for the normal fluid component (velocity v n ) and no- 
penetration for the superfluid component (veloci ty v, 9 ), ignoring; the sma ll te nsion force to 



reduce the order of the equation for v s by one [see iHenderson et al.l (119951 ) and lPeralta et al. 



(120081 ) for a discussion] . The mutual friction force is taken to h ave the anisotropic H all-Vinen 



form [oc tb s x <jj s x v ns , with v 
B = 10~ 2 , and B' 



and u = V x v s jHall fc Vinenlll956al lb[)]. with 
10- 4 (lAndersson et alil2007h . We take p n /p = 0.01 and p s /p = 0.99 
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for the normal an d superfluid mass d ensity fractions respectively, where p 
the total density (IPeralta et al.ll2005l ). 



p n + p s denotes 



3. Nonaxisymmetric instabilities 

The Stewartson layer is disrupted when Ro exceeds a thresh old Ro r {E). which de- 
creases as E increas e s, exc iting a Kelvin-Helmholtz-type instability. IHollerbachl (120031 ) and 
Schaeffer &: Cardinl (120051 ) investigated thoroughly the most unstable modes of a viscous 
fluid Stewartson layer. They discovered nonaxisymmetric instabilities for azimuthal modes 
1 < m < 119, with < Ro c < 0.6 and 10 -10 < E < 10~ 3 5 . The azimuthal mode number of 



the most unstable mode increases with decreasing E. IHollerbachl (120031 ) discovered empiri- 
cally the scaling Ro c ~ 0.6E om for Ro > 0. For Ro < 0, the most unstable mode is almost 
always m = 1, wi th m = 2 being excited in the range 10~ ' 25 £ Ro c ;$ 0.1. In this regime, 
Hollerbachl (120031 ) foun d Ro r ~ 0.8E OA 5 . The asymmetry with the sign of Ro is still not 



understood physically (IHollerbachl 120031 1 . 



Here, we generalize the numerical calculation of Ro c to a superfluid Stewartson layer. 
We follow a two-stage recipe. First, for fixed Ro = 10~ 4 , we calculate axisymmetric HVBK 
basic states for E = 1 x 10~ 3 , 1 x 10~ 4 , 5 x 10~ 5 , 2 x 10" 5 , and 1 x 10~ 6 . Second, we 
linearize the HVBK equations around the base state and test the stability of a given m/fl 
perturbation (typical amplitude ~ 1 %) by increasing Ro until the mode grows exponentially. 
We obtain the scaling 



Ro r « 4.1E 



0.40 



for 10~ 6 < E < 10~ 3 and 0.02 < Ro c < 0.33. We concentrate our efforts on m = 6, the most 
unstable mode at E w 10~ 3 for viscous fluids (and also for the superfluid). It is important to 
note that m = 6 is not the most unstable mode at smaller Ekman numbers (e.g. E = 1 x 10~ 6 , 
where m = 10 is more unstable), but the critical Rossby number is unaffected (Ro c = 0.02). 
In $3j w e extrapolate the scaling (ED) to the regim e E £ 10~ 12 , Ro £ 10~ 4 , relevant to radio 



pulsars (ILyne et al.ll2000l ; iMelatos fc Peraltall2007l ). where numerical simulations are too hard 
to do. A more thorough parameter survey will be presented in a future paper. 

Nona xisymmetric instabilitie s can decrease the shear inside a viscous Stewartson layer 
by 30 % (IHollerbach et al.l 120041 ) . We observe a similar but less pronounced effect in the 
superfluid, e.g. when the mode m = 6 is excited at Ro c = 0.3, with E = 10~ 31 . Figure 
[2] describes the topology of the flow in the layer, before (Figure [2^, Ro = 0.2) and after 
(Figure [2b, Ro = 0.3) the mode m = 6 is excited. The discriminant Da = Q\ + 27/^/4, 
with Ra = det(Aij), Qa = (Af, — AijAji)/2, and A±j = dvi/dxj, distingui shes between 



regions that are focal (Da > 0, blue) and strain-dominated (Da < 0, green) (jChong et al. 
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19901 ). The maximum shear changes from doo/ds — 1.1 at t — to du/ds = 1.0 at t = 400 



(cf. viscous fluid, where the observed change is from du/ds = 1.3 to 0.8). The hexagonal 
flow structure also boosts the torque one must exert on the inner sphere to maintain the 
shear, although the increase is less dramatic (~ 10 %). 



4. Pulsar glitches 



As E and Ro control the stability of the Stewartson layer, it is interesting to test whether 
the amplitude and rate of incidence of rotational irregularities like glitches are related to these 
two dimensionless quantities. In order to calculate E, we need to know how v n depends on 
the density p and temperat ure T in the outer core. U sing the neutron-neutron scattering 
viscosity formula derived by lCutler fc Lindbloml (119871 ). we find 



E = 2.6 x 10 



-12 



T 



10 s K 



10 2 rads - 



(2) 



2.8 x 10 12 g cm 3 . The core temperatur e T is related to the s urfac e temperature 



with p 

T s , e.g., via the two- zone heat-blanket model of iGudmundsson et al.l (119821 ). which gives 
T/10 8 K = 1.29(T S /10 6 K) 18 . We estimate T s from the characteris tic age r c = n 2 /2\tt 2 \, 
combined with theoretical cooling curves for r c < 10 6 yr ( jPagdll998l ) and standard neutrino 
cooli ng; similar E distributions are obtained with nonstandard cooling ( IMelatos &: Peralta 
2007h . 



In Figure [31 we plot a s points the m a ximum Ro in 55 glitching pulsars with r c < 10 6 
yr, taken from Table 1 in IMelatos et al.l (120081 ) . assuming conservatively that Ro is less 
than the fractional frequency jump of the largest glitch. For each object, E is estimated 
by the method in the previous paragraph. We also plot two Ro c (E) curves: the HVBK 
scaling c omputed in jj2] (Ro r ~ 4.1 E 10 - 40 , solid curve) and a scaling extrapolated from the 
study by ISchaeffer fc Cardinl J2005h for Ro c < 10" 2 and E < 10" 5 (Ro c w g.4^ - 57 , dashed 



curve). The viscous fluid scaling lies in the middle of the cluster of points, while the HVBK 
superfluid scaling lies above all the points. This suggests two possible conclusions: (i) a 
glitching pulsar must have Ro < Ro c (E), otherwise Stewartson layer instabilities would 
erase the shear required for the glitch phenomenon to occur; and (ii) it is important to 
include the HVBK superfluid dynamics to ensure that all points in Figure [3] lie below the 
Ro c (E) curve, given that no discernible difference is observed in the glitch behaviour of 
objects above and below the viscous fluid scaling. 



In the superfluid Stewartson layer, there is an extra ingredient that influences the insta- 
bility curve: the mutual friction between the normal and superfluid component, controlled 
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by the dimensionless parameters B and B' . We have only considered the weak coupling 
regime i n this investigation (B = 10 -2 , B' = 10~ 4 ). In the strong coupling regime [B = 0, 
B' = 1 (lAndersson et al.ll20071 )] the instabilities are likely to be quite different. Preliminary 
results for axisymmetric steady states show that the Stewartson layer is ~ 10 % thicker. A 
more detailed investigation of the effect of mutual friction and entrainment on Ro c (E) will 
be presented in a future paper. 



5. Discussion 



The results in §2] and §0 demonstrate that Stewartson layers develop in rotating HVBK 
superfluids, like the 1 So-paired neutron superfluid in the outer core of a neutron star. We 
present the first numerical simulation of a HVBK Stewartson layer, for Ekman and Rossby 
numbers in the ranges 10~ 6 < E < 10~ 3 and 10 -4 < Ro < 10~ 2 . Superfluid Stewartson layers 
are unstable to nonaxisymmetric perturbations. Transitions between unstable Stewartson 
modes in the superfluid are different to those found in a viscous fluid, and the critical Rossby 
number Ro c (E) is ~ 10 times higher in the above parameter regime. In glitching pulsars, 
one finds Ro < Ro c (E) in the 55 pulsars for which Ro and E can be estimated reliably from 
observations. One possible interpretation of these data is that the Stewartson layer remains 
stable in these objects, allowing rotational shear to build up (as required for glitches) without 



trigg ering disruption of the Stewartson layer [which would nullify the shear (IHollerbach et al. 



2004 ) and hence shut down the glitching behaviour]. The threshold Ro c (E) can be compared 
against the upper limit derived independently from the gravitational- wav e spin down caused 
by K olmogorov-like superfluid turbulence excited in the stellar interior (IMelatos fc Peralta 
2009h . 



The conclusions drawn from the data in Figure [3] extend and partially clarify the sur- 
prising empir ical finding, that the E d istribution is markedly different between glitchers and 
nonglitchers (IMelatos fc Peraltal 120071 ). It seems strange that stars with E ~ 10~ 10 (few 
glitches) and E ~ 10~ 12 (many glitches) behave so differently, since in both regimes Kol- 
mogorov t urbulence must be fully developed and scale-free. However, from Figure [3] [and 
Figure 4 in IHollerbach! (120031 )] one sees t hat small differences in E and Ro lead to very differ- 
ent flow states and stability properties (lHollerbachll2003l ; IHollerbach et al.ll2004l ). Moreover, 
if theoretical estimates of effective viscosity and hence E are too low by a factor ~ 10 5 , 
due to turbulent Reynolds stresses (IMelatos fc Peraltal 120071 ). glitching pulsars lie in a range 
where (i) the most interesting flow transitions occur before the flow becomes turbulent (i.e. 
at 10~ 4 E ^ 10 -2 for 1 < m < 6 ), and (ii) numerical simulations are computationally 
tractable. 
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The results of this paper do not prove that Stewartson flow transitions control glitch 
behaviour. We merely find empirically that all observed glitchers lie on the stable side of the 
Ro c (E) threshold for nonaxisymmetric instabilities of a Stewartson layer in a HVBK super- 



fluid. Recently, a hydrodynamic trigger for pulsar glitches was proposed by lGlampedakis Andersson 



( 120091 ) . associated with r- modes excited by differential rotation. It would be interesting to 
see how meridionally circulating Stewartson base states modify these calcul ations, especially 
in the strong pinning scenario (for which B' = 1, B = 0). 'Price et al.l (120081 ) found empirical 
evidence of departures from solid-body rotation in a radio pulsar. 

We do not consider stratification or magnetic fields in this investigation. In viscous 
fluids, strong stratifi cation supresses the Stewartson layer, leaving a Taylor column parallel 
to the rotation axis (IHollerbachl l2009f). M agnetic fields widen the Stewartson layer, merging 
it with the interior flow (IHollerbachl 1 19941 ). This difficult physics deserves investigation. 
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(a) 




Fig. 1. — (a) Angular velocity u as a function of cylindrical radius rsin#, at z — 1.2, for 
the normal fluid component in superfluid spherical Couette flow, for E — 1 x 10~ 3 (upper 
curve), 1 x 10~ 4 (middle curve), and 2 x 10~ 5 (lower curve), (b) Meridional streamlines of 
the normal component for E — 1 x 10~ 3 (left), 1 x 10~ 4 (center), and 2 x 10~ 5 (right). Ro 
equals 10~ 4 in all plots. 




z 



Fig. 2. — Isosurfaces of the velocity gradient discriminant (see text): Da = lO -3 ^ (blue) 
and — 10~ 3 ^2 (green) for E = 10~ 31 and the transition Ro = 0.2 — > 0.3, where m = 6 is the 
most unstable mode, (a) t — (before transition); (b) t = 400 (after transition). 
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Fig. 3. — Critical Rossby number Ro c versus Ekman number E for 55 glitchers with charac- 
teristic ages r c < 10 6 yr. The points indicate upper limits on Ro c derived from the largest 
observed glitch in each object. E is calculated using standard neutrino cooling. The solid 
curve is a fit to the HVBK superfluid simulations i n jj2] and jj3l The dashed curve is an 
extrapolated fit to the viscous fluid calculations from ISchaeffer fc Cardinl (120051 ). 



